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EFFECTS  OF  CAVITATION  ON  UNDERWATER  SHOCK  LOADING 
-  PLANE  PROBLEM,  FINAL  REPORT 

1.  Introduction 

Earlier  phases  of  this  investigation  were  described  in 
a  series  of  reports.  Problems  of  dependent  variable  selection 
were  considered  in  Ref.  1.  The  viability  and  utility  of  a 
displacement  potential  were  demonstrated  by  solving  a  repre¬ 
sentative  one-dimensional  problem  and  obtaining  results  in 
good  agreement  with  those  obtained  earlier  by  Bleich  and 
Sandler  and  reported  in  Ref.  2.  An  unsuccessful  effort  to 
use  ADINA  for  applications  with  axisymmetric  geometry  was 
described  in  Ref.  3. 

Development  of  a  finite  element  FORTRAN  IV  program 
(DPLPOT)  to  trace,  for  plane  regions,  the  propagation  of 
underwater  shock,  including  effects  of  response-induced 
cavitation,  was  described  in  Ref.  4.  A  mesh  generator  (MGNEW) 
was  also  developed  for  these  applications.  Reference  5 
describes  the  construction  of  a  computer  program  (STRUK5)  to 
determine  the  transient  response  of  a  submerged  cylindric 
structure.  Early  results  obtained  with  the  concurrent  use  of 
STRUK5,  or  its  successor,  STRUK6 ,  and  the  fluid  program  DPLPOT 
were  reported  in  Refs.  6  and  7. 

The  present  report  summarizes  results  for  the  plane 
problem.  Data  presented  are  representative  of  findings  from 
over  500  separate  computer  runs.  For  the  range  of  parameters 
considered,  conclusions  are  given  concerning  the  effects  of 
cavitation  on  the  severity  of  structural  loading. 


2.  Scope  of  Investigation 
2 . 1  Fluid  Parameters 

Properties  of  the  fluid  are  representative  of  seawater 
with  density  *  1024  kg/m3  and  acoustic  velocity  c  *  1500  m/s. 
The  fluid  region  is  rectangular  (see  Fig.  1) .  Hydrostatic 
pressure  P^  is  uniform  because  a  gravity  induced  gradient 
would  violate  the  assumed  symmetry  relative  to  y  *  0. 

Incident  shock  waves  have  plane  fronts  (x  *  const.) 
with  step  rise  of  height  P  ,  followed  by  exponential  decay 
with  characteristic  length  L.  To  simulate  the  effect  of 
surface  cutoff,  a  limited  number  of  runs  were  made  using 
box  waves  of  length  L. 

The  piecewise  linearity  of  the  fluid  constitutive  law 
(Ref.  1)  reduces  the  number  of  different  combinations  of  P^ , 
Pg,  and  L  that  need  to  be  investigated.  Specifically,  if  L 
and  Pg/Pfj  are  fixed,  response  variables  such  as  shell  bending 
moments,  stresses,  and  accelerations  vary  in  direct  propor¬ 
tion  to  the  shock  pressure  Fs. 

The  shell  radius  used  throughout  is  R  =  5  m.  The  decay 
lengths  of  incident  shock  waves  were  7.5  m,  15  m,  and  30  m, 
corresponding  to  decay  times  of  5  ms,  10  ms,  and  20  ms.  As 
reported  previously,  when  cavitation  occurs  it  begins  at 
least  one  decay  length  in  front  of  the  structure  (x  >  R  +  L 
in  Fig.  1) .  Since  the  boundary  conditions  at  entry,  top,  and 
exit  faces  are  valid  only  if  the  adjacent  fluid  is  uncavi- 
tated,  the  required  region  size  is  problem  dependent.  Except 
in  very  short  runs  used  to  look  at  fine  details  of  the  early 
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behavior,  the  exit  face  was  located  at  x  »  -13.5  m.  Position 
of  the  top  face  ranged  from  y  *  22  m  to  y  *  72  m.  For  the 
entry  face,  the  normal  range  was  from  x  =  36  m  to  x  =  108  m. 

A  useful  program  feature  that  allows  determination  of 
cavitation  effects  by  direct  comparison  is  an  "ON-OFF" 
switch  for  cavitation.  In  the  "OFF"  mode,  the  fluid  pressure 
is  the  same  linear  function  of  condensation  for  negative 
pressure  as  it  is  for  positive  pressure  and  no  cavitation 
occurs . 

2. 2  Structure  Parameters 

As  noted  above,  a  shell  radius  R  »  5  m  was  used  through¬ 
out.  The  sandwich  shell,  made  from  steel  with  modulus  of 
elasticity  E  *  210  GPa  and  density  pg  *  7830  kg/m3,  consisted 
of  inner  and  outer  shells,  each  25  mm  thick,  joined  by  a  mass¬ 
less  core  291  mm  deep.*  The  core  is  assumed  to  make  no  contri¬ 
bution  to  circumferential  or  flexural  rigidity,  but  to  provide 
a  rigid  shear  connection  between  inner  and  outer  shells. 

The  addition  of  a  single  internal  mass  ("contents"  mass), 
elastically  coupled  to  the  shell  rigid  body  mode  (n  *  1) ,  was 
described  in  Ref.  6.  A  subsequent  refinement  has  added  viscous 
damping  in  parallel  with  the  elastic  suspension  of  the  contents 
mass.  Because  the  contents  mass  has  its  magnitude  adjusted 
to  make  the  assembly  neutrally  buoyant,  the  only  additional 
parameters  introduced  by  its  inclusion  are  the  (fixed  base) 


*This  unlikely  core  depth  corresponds  to  (r  /R) 2  *  .001, 
where  r^  is  the  section  radius  of  gyration? 
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frequency  f  and  critical  damping  fraction  z  of  the  suspension. 
For  most  of  the  cases  investigated  the  choices  f  =  10  Hz  and 
5  =  .10  have  been  made. 

2.3  Spatial  and  Temporal  Subdivision 

Appropriate  choices  of  mesh  subdivision  and  the  step  size 
used  in  time  integration  have  been  the  subject  of  extensive 
investigation.  It  is  clear  that  finer  grids  are  capable  of 
producing  improved  accuracy.  If  both  structure  and  fluid  were 
governed  by  elliptic  partial  differential  equations,  the  ideal 
compromise  would  provide  a  fine  grid  at  the  structure  and  in 
the  adjacent  fluid.  In  the  remainder  of  the  fluid,  the  element 
size  would  increase  with  distance  from  the  structure.  Because 
the  fluid  behavior  is  governed  by  a  hyperbolic  partial  differ¬ 
ential  equation,  such  a  graded  mesh  is  not  satisfactory. 

The  dilemma  posed  is  as  follows.  Satisfactory  simulation 
of  wave  propagation  in  the  fluid  requires  a  substantially  uni¬ 
form  subdivision  throughout.  The  choice  made  for  the  fluid 
determines  the  number  of  interface  fluid  nodes  and  compati¬ 
bility  requirements  at  the  interface  dictate  an  equal  number 
of  structural  nodes.*  Thus,  the  structural  subdivision  is 
fixed  by  the  choice  made  for  the  fluid.  Best  accuracy  for 
time  integration  in  the  fluid  requires  the  largest  time  step 
compatible  with  stability,  whereas  the  accuracy  attained  for 
the  structure  improves  as  the  time  step  size  is  decreased. 


*This  constraint  results  from  the  kind  of  structural  model 
used  here  and  is  not  a  necessary  condition  for  all  fluid- 
structure  interactions. 


In  the  face  of  these  conflicting  requirements,  the 
compromise  adopted  was  based  on  the  following  choices. 
For  full  field  problems,  the  number  of  fluid  nodes  and 
the  number  of  fluid  elements  was  approximately  3000. 

The  time  step  was  chosen  to  be  approximately  one-fifth 
the  period  of  the  highest  frequency  structural  mode. 
Mesh  data  and  time  steps  are  summarized  in  Table  1  for 
the  three  decay  lengths  investigated. 


Table  1.  Mesh 

Data 

Decay 

length , 

L 

7.5m 

15  m 

30  m 

Number  of  fluid  nodes 

3558 

3088 

3004 

Number  of  fluid  elements 

3392 

2950 

2880 

Element  side,  largest 

.69  m 

1.11  m 

1.84  m 

Element  side,  smallest 

.49  m 

.79  m 

1.31  m 

Structural  degrees  of  freedom 

64 

40 

24 

Time  step 

.031  ms 

.05  ms 

.10  ms 

8 


3.  Response  as  a  Function  of  Time 


3 . 1  Parameters  for  Chosen  Cases 

In  the  following  sections,  time  histories  for  shell 
bending  stress,  hoop  compressive  stress,  and  contents  accel¬ 
eration  are  presented  for  three  selected  encounters.  Parame¬ 
ters  used,  together  with  Figure  numbers  for  the  relevant 
curves,  are  given  in  Table  2. 

Table  2.  Parameters  for  Time  Histories 


L 

?h 

PS 

9 

Figure  Numbers 

7.5  m 

.4  MPa 

12  MPa 

0° 

2, 

5,  8 

15  m 

.4  MPa 

12  MPa 

180° 

3, 

6,  9 

30  m 

.  2  MPa 

12  MPa 

0° 

4, 

7,  10 

Time  is  measured  from  the  instant  of  first  contact  between 
shock  front  and  shell.  This  occurs  when  the  shock  front 
is  at  x  =  R  =  5  m. 

3 . 2  Shell  Bending  Stress 

At  the  location  of  first  contact  (9  =0,  Fig.  1) ,  the 
bending  stress  rises  rapidly,  usually  attaining  its  first 
maximum  within  5  ms.  A  representative  set  of  curves  showing 
bending  stress  as  a  function  of  time  is  given  in  Fig.  2. 
Positive  bending  stress  corresponds  to  compressive  stress  on 
the  external  surface  of  the  hull.  Results  are  for  decay 
length  L=7.5mat9=0.  The  curve  labeled  ON  reveals  that 
cavitation  causes  a  higher  maximum  bending  stress  at  T  *  68  ms 
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than  at  T  =  1  ms.  The  curve  labeled  OFF  (cavitation  suppressed) 
shows  the  same  early  peak,  but  no  comparable  late-time  maxi¬ 
mum.  A  parallel  set  of  results  for  9  =  180°  and  L  =  15  m  is 
given  in  Fig.  3.  A  third  set,  for  9=0  and  L  =  30  m,  is 
given  in  Fig.  4. 

In  order  to  assess  the  effects  of  cavitation  on  stress, 
it  is  necessary  to  know  where  and  when  cavities  appear. 

In  all  cases  based  on  realistic  parameters,  the  first  cavity 
appears  at  a  distance  of  about  one  decay  length  L  in  front 
of  the  hull  after  at  least  one  decay  time  (T  =  L/c) * . 

This  cavity  and  the  effects  of  its  closure  do  not  percep¬ 
tibly  influence  shell  bending  stresses.  At  a  later  time, 
there  may  be  cavitation  at  and  near  the  hull.  The  late 
time  bending  stress  extrema  in  Figs.  2,  3,  and  4  all  coin¬ 
cide  approximately  with  closures  of  such  hull  cavities. 

3 . 3  Hoop  Compressive  Stress 

Superposed  on  the  circumferential  compressive  stress 
resulting  from  external  hydrostatic  pressure  are  very  high 
transient  stresses  induced  by  the  shock  loading.  Represent¬ 
ing  compressive  stresses  as  positive,  dynamic  hoop  stresses 
as  a  function  of  time  are  given  in  Figs.  5,  6,  and  7  for 
the  same  cases  as  in  Figs.  2,  3,  and  4. 

In  all  cases,  it  will  be  observed  that  the  first  direct 
stress  extremum  occurs  between  T  =  7  ms  and  T  =  11  ms.  This 


*For  the  box  wave,  cavitation  first  occurs  about  one-half 
box  length  in  front  of  the  hull  and  at  an  earlier  time. 


10 


extremum  is  not  affected  by  cavitation.  At  late  times, 
the  effects  of  cavitation  are  discernible,  but  the  corres¬ 
ponding  extrema  are  no  greater  than  15%  of  the  initial  one. 
3.4  Contents  Acceleration 

The  acceleration  of  the  contents  mass  is  an  indication 
of  the  severity  of  the  shock  motion  transmitted  to  equipment 
Time  history  of  the  contents  acceleration  is  presented  in 
Figs.  8,  9,  and  10  for  the  three  cases  previously  considered 
.e  first  extremum  occurs  before  T  =  12  ms.  The  second  ex¬ 
tremum  is  usually  of  approximately  equal  magnitude,  but 
opposite  sign.  The  effects  of  cavitation  upon  the  largest 
(absolute)  accelerations  are  insignificant. 
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4.  Structural  Importance  of  Cavitation 


4 . 1  Bending  Stress 

The  hull  bending  stresses  shown  in  Figs.  2,  3,  and  4 
all  show  late  time  extrema  associated  with  the  closures  of 
cavities  at  the  hull.  In  Fig.  3,  this  peak  is  larger  and 
of  greater  duration  than  the  early  one.  It  is  believed 
that  the  hull  cavitation  encountered  is  induced  by  the 
large  amplitude  shell  motion  that  occurs  in  phase  opposi¬ 
tion  to  the  response  motion  of  the  contents  mass.  This 
behavior  was  illustrated  in  Figs.  1  and  2  of  Ref.  6  and  is 
shown  again  here  as  Fig.  11.  A  more  realistic  and  elabor¬ 
ate  model  of  the  contents  would  introduce  many  additional 
degrees  of  freedom  with  a  diversity  of  individual  frequen¬ 
cies.  The  accompanying  shell  motion  might  reasonably  be 
expected  to  have  smaller  amplitude  and  be  less  regular. 

An  indication  of  the  dependence  of  the  hull  cavitation 
effects  on  the  contents  model  is  shown  in  Fig.  12.  Here, 
the  solid  curve  is  the  same  as  in  Fig.  3  and  the  dotted 
curve  shows  what  happens  to  shell  bending  stress  when  the 
contents  suspension  frequency  is  reduced  from  10  Hz  to  5  Hz. 
This  change  eliminates  the  hull  cavitation  and,  thus,  the 
late  time  stress  peak.  There  is  a  striking  similarity  be¬ 
tween  the  dotted  curve  of  Fig.  12  and  that  of  Fig.  3.  In 
the  latter  case,  cavitation  was  artificially  suppressed. 

For  the  parameters  of  Fig.  2,  a  change  of  suspension 
frequency  from  10  Hz  to  5  Hz  will  again  suppress  hull  cavita¬ 
tion  and  eliminate  the  late  time  stress  peak.  The  effect 
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is  shown  in  Fig.  13. 


In  consideration  of  these  results,  it  is  concluded  that 
the  late  time  bending  stress  peaks  of  Figs.  2,  3,  and  4  are 
spurious  effects  of  the  oversimplified  contents  model.  It 
follows  that  cavitation  does  not  increase  the  extreme 
values  of  bending  stress. 

4 . 2  Hoop  Compressive  Stress 

The  results  presented  in  Figs.  5,  6,  and  7  show  that 
early  time  peaks  in  dynamic  hoop  compressive  stress  are 
unaffected  by  cavitation.  Although  cavitation  does  seem  to 
make  late  time  peaks  higher,  the  latter  are  so  small  in 
comparison -with  the  early  peaks  that  the  cavitation  effects 
do  not  affect  structural  survival. 

Attention  is  directed  to  the  very  large  magnitudes  of 
peak  values  of  hoop  compressive  stress  compared  with  those 
associated  with  shell  bending.  It  is  clear  that  the  largest 
resultant  stresses  occur  at  early  times  and  that  late  time 
bending  stress  peaks  will  produce  much  smaller  resultants 
than  at  early  times.  This  lends  additional  support  to  the 
proposition  that  cavitation  does  not  increase  the  likeli¬ 
hood  of  hull  collapse. 

4. 3  Internal  Equipment 

The  contents  mass  of  the  structural  model  may  be  con¬ 
sidered  to  represent  a  lumping  of  internal  structure  and 
equipment.  It  is  a  very  crude  model,  but  the  contents 
accelerations  presented  in  Figs.  8,  9,  and  10  show  that 


extreme  values  are  not  significantly  affected  by  cavitation. 
Although  a  more  elaborate  model  would  be  required  to  predict 
with  confidence  the  accelerations  of  particular  equipment 
items,  the  conclusion  that  cavitation  does  not  significantly 
affect  peak  structural  accelerations  is  unlikely  to  be 
reversed  by  using  a  more  complex  model. 

4 . 4  Effects  of  Surface  Cutoff 

At  moderate  submergence  depths  with  large  horizontal 
standoff,  the  incident  shock  pulse  is  abruptly  truncated 
by  a  negative  pressure  wave  reflected  from  the  free  surface. 
As  was  mentioned  above,  a  limited  number  of  runs  were  made 
using  an  incident  box  wave  to  simulate  this  effect.  Details 
of  response  to  this  altered  wave-form  are  omitted  here  be¬ 
cause  the  behaviors  found  are  qualitatively  identical  with 
those  reported  in  Section  3 . 


Conclusions 


5 . 

For  the  response-induced  cavitation  that  has  been  the 
subject  of  this  investigation,  none  of  the  results  obtained 
indicates  that  cavitation  increases  the  likelihood  of  struc¬ 
tural  or  equipment  failure.  There  is  a  remaining  possibility 
that  the  large  deflections  accompanying  inelastic  shell 
response  might  induce  close-in  cavitation  whose  collapse  would 
produce  higher  extreme  structural  loadings.  Study  of  in¬ 
elastic  response  and  the  accompanying  geometric  nonlinearity 
is  beyond  the  scope  of  this  investigation. 


6.  Reconmendat ion 


The  investigation  which  is  completed  with  this  report 
has  dealt  with  response-induced  cavitation.  The  parameters 
have  been  chosen  to  represent  attacks  on  submerged  submarines 
and  bulk  cavitation  induced  by  the  free  surface  has  been  ex¬ 
cluded. 

There  is  abundant  evidence,  both  experimental  and 
analytic,  that  bulk  cavitation  and  the  subsequent  cavity  col¬ 
lapse  may  produce,  for  surface  ships,  significant  increases 
in  the  severity  of  shock  loadings.  It  is  recommended  that 
the  finite  element  technique  based  on  a  potential  function 
be  extended  to  the  important  class  of  problems  featuring 
water  borne  shock  loads  in  the  presence  of  bulk  cavitation. 
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